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Abstract. This programmatic paper lays down the possibility to reconcile the necessity to resum many-body correla¬ 
tions into the energy kernel with the fact that safe multi-reference energy density functional (EDF) calculations cannot 
be achieved whenever the Pauli principle is not strictly enforced, as is for example the case when many-body correla¬ 
tions are parametrized under the form of empirical density dependencies. Our proposal is to exploit a newly developed 
ab initio many-body formalism to guide the construction of safe, explicitly correlated and systematically improvable 
parametrizations of the off-diagonal energy and norm kernels that lie at the heart of the nuclear EDF method. The many- 
body formalism of interest relies on the concepts of symmetry breaking and restoration that have made the fortune of 
the nuclear EDF method and is, as such, amenable to this guidance. After elaborating on our proposal, we briefly outline 
the project we plan to execute in the years to come. 

PACS. 21.60.Jz Nuclear DFT and extensions - 21.30.Cb Nuclear forces in vacuum - 21.30.Fe Forces in hadronic 
systems and effective interactions - 21.60.De Ab initio methods 


1 Introduction 


Nuclear ab initio methods combine state-of-the-art models of 
elementary interactions with controlled many-body expansion 
techniques. Confronting their results with experimental data, 
ab initio calculations provide a test of our understanding of the 
strong interaction between nucleons, at least as long as the er¬ 
ror associated with the solving of the many-body Schrodinger 
equation is smaller than the residual distance to the data. In 
that sense, such methods are not meant to account for data at all 
cost but rather use the discrepancy with them to gauge the qual¬ 
ity of the input. When both nuclear interaction models and the 
solving of the many-body Schrodinger equation reach a good 
enough accuracy, ab initio methods can provide controlled ex¬ 
trapolations to experimentally unknown regions. 

Traditionally limited to the lightest nuclei [l||2][3|, ab initio 
many-body methods have been extended tremendously over the 
last ten years to medium-mass closed-shell nuclei 05][6l7f8j 
|9} and to those displaying an open-shell character B1011111112] 
IT5IT51 . Based on this extended reach, ab initio many-body cal¬ 
culations are providing an even more stringent test of nuclear 
interactions IU51I131116M built within the frame of, e.g., chiral 
effective-field theory [17)181 than they have done so far. 


The nuclear energy density functional (EDF) method pro¬ 
vides a powerful quantum mechanical tool that can be applied 
to all bound atomic nuclei lfl9l . irrespective of their mass and 
isospin, thanks to a low computational cost. The EDF method 
is thus characterized by an extended reach and aims, from the 
outset, at accounting for empirical phenomena with the highest 
possible precision, at least in the vicinity of the region where 
experimental data are available to adjust the parameters enter¬ 
ing the EDF kernels. It is, however, an effective method that 
does not provide (direct) information about interactions be¬ 
tween elementary nucleons. Most importantly, the empirical 
formulation of currently available parametrizations of the EDF 
kernels, i.e. the fact that they are not rooted into sound many- 
body methods, makes EDF predictions away from known data 
unreliable if not plagued with spuriousities leading to criti¬ 
cal pathologies li20l . Those pathologies are the consequences 
of a violation of the Pauli exclusion principle that eventually 
contaminates state-of-the-art multi-reference EDF calculations 
with nonphysical contributions to the energy l2Tll22ll23l . The 
violation of the Pauli principle in the EDF kernels relates itself 
to the traditional way of parametrizing many-body correlations 
under the form of density dependencies or to a relaxation of 
specific interrelations between the coefficients of the various 
terms at play. Solutions to better formulate the multi-reference 
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EDF method and in particular the restoration of symmetries 
constitute a top priority today I124II25II26) . 

It is the goal of this programmatic document to propose 
one possible way to overcome these limitations by rooting the 
formulation of the nuclear EDF method into sound and appro¬ 
priate many-body techniques. As sketched in Fig. Q] there are 
typically two ways how ab initio many-body methods can be 
used to improve on the current status of EDF calculations 

1. When ab initio methods based on realistic inter-nucleon in¬ 
teractions are mature enough, their predictions for exper¬ 
imentally unknown nuclei can be used as pseudo-data to 
better constrain parameters entering the functional form of 
the nuclear EDF kernels. As of today, however, ab initio 
calculations of mid-mass nuclei, which constitute the natu¬ 
ral overlap region with the EDF method, have not reached 
such a maturity yet. 

2. The expansion of the Schrodinger equation at play in a 
given many-body method can be used as a mathematical 
guidance to build sound parametrizations of the EDF ker¬ 
nels that incorporate much needed correlations while avoid¬ 
ing unwanted pathologies. For this rationale to be operative, 
the many-body technique of interest must build on the same 
key concepts as those underlying the EDF method, i.e. the 
breaking and the restoration of symmetries. It is only in 
the last five years that ab initio many-body methods com¬ 
bining self-consistent Green’s function (SCGF) ll llll3l or 
coupled-cluster (CC) 110112711141 with the concept of sym¬ 
metry breaking have been implemented to tackle open-shell 
nuclei. It is even more recently that methods adding the ex¬ 
act restoration of the broken symmetry have been formu¬ 
lated 1281 . 

While the first point must be postponed to several years in the 
future, the present work elaborates on the second. The proposal 
is fundamentally based on the novel many-body formalism, or 
more specifically on its simplified many-body perturbation the¬ 
ory (MBPT) version, proposed in Ref. (28). As will be de¬ 
scribed below, this formalism consistently builds on the suc¬ 
cessive breaking and restoration of symmetries and is, as such, 
the first and only full-fledged many-body technique that can 
provide a constructive approach to off-diagonal EDF kernels. 

The paper is organized as follows. Section [2] provides a 
brief account of the nuclear EDF formalism as it stands to¬ 
day. Basic ingredients, key concepts and current limitations are 
underlined to facilitate in Sec. [3] the introduction of the novel 
many-body method used as a guidance to build the new fam¬ 
ily of EDF parametrizations. The extended EDF scheme that 
emerges from this proposal is outlined in Sec. [4] Conclusions 
and perspectives are given in Sec. [3] 

2 Modem nuclear EDF formalism 

We provide a brief introduction to the nuclear EDF formal¬ 
ism based on the more complete account given in Ref. ||26) . 
As such, we do not aim at reviewing the status of the field, at 
covering all possible ramifications of the approach or at pre¬ 
senting recent achievements and applications. For standard re¬ 
views that cover the connection to empirical data, we refer the 
reader to, e.g.. Refs. Ill9ll29l . 



Fig. 1. (Color online) Potential connection and cross fertilization be¬ 
tween EDF and ab initio many-body methods. 

2.1 Reference states and off-diagonal EDF kernels 

The ingredients lying at the heart of the EDF method are the 
so-called off-diagonal norm N(g';g ) and energy ll(g'\g) ker¬ 
nels ll26l . The norm kernel is traditionally defined as the plain 
overlap 

N(g'-g) = (<P(g t )\<P(g)), ( 1 ) 

between two many-body states of Bogoliubov type. The latter 

designates normalizecQ product states of the form 

!*(*)> = *m w l 0 >. (2) 

n 

where quasi-particle creation and annihilation operators satis¬ 
fying {j3relate to particle operators {c„;c a } 
associated with a basis of the one-body Hilbert space through 

the so-called Bogoliubov transformation 

Pjf ) =Z(ujft*c a +vj*>*cl l ), (3a) 

a ' 7 

/3^ )+ = E (Va g Jc a + Ui g >c a ) . (3b) 

a v 7 

Matrices U^ and V^ make up the unitary Bogoliubov trans¬ 
formation If30l . 

The collective index g = \g\e’ n labeling the many-body 
states gathers a set of order parameters characterizing the po¬ 
tential breaking of symmetries of the underlying nuclear Hamil- 
toniarQ. The norm |g| of the order parameter tracks the extent 

1 This corresponds to using the convention N(g\g) = 1. 

2 In the most general setting, the label g may also incorporate non¬ 
collective quantum numbers characterizing a set of quasi-particle ex¬ 
citations 03332). 
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to which |<f>(g)) breaks the symmetry, i.e. its ’’deformation”, 
whereas the phase £2 = Arg(g) characterizes the orientation of 
the deformed body with respect to the chosen reference framtQ- 
For nuclei, the main symmetries to be considered are (i) SU (2) 
associated with rotational invariance in real space and the con¬ 
servation of angular momentum and (ii) t/(l) associated with 
rotational invariance in gauge space and the conservation of 
particle number. The associated compact Lie group is generi- 


cally defined as Sf = {R(f2);f2 £ Dy} and the 
tween states differing by their angles writes as 

connection be- 

\<P(\g\,£2))=R(£2)\cP(\glO)). 

(4) 

In general, the off-diagonal energy kernel 


H{g'\g)=h{g'\g)N{g'\g), 

(5) 

invokes a functional 


h{g'-,g)=h[p s ' g , K*'*,K**'*] 

( 6 ) 


of the off-diagonal normal and anomalous one-body density 
matrices constructed from the two product states involved through 


g' g (^(sOkyCil^g)} 

PiJ ~ W)l*(*)> ’ 

(7a) 

g'g (^fe')|C;C,|^(g)) 
ij ~ W)\<P(g)) ’ 

(7b) 

g gU {^(s')\c f ic]\<P(g)) 

K U ~ (0(g')\0(g)) ’ 

(7c) 


The two reference states implicated |<£>(g)) and ^(g')) are a 
priori different and provide the kernels with their off-diagonal 
character. 


2.2 SU(2) symmetry as an example 

For presentation purposed we limit the discussion throughout 
the rest of the paper to the SU (2) = {R(f2),f2 £ D SU ( 2 )} group, 
where £2 = (a,/3,y) actually embodies the three Euler angles 
whose domains of definition are 


D SU ( 2 ) — D a x Dp x Dy = [0,4zr] x [0, tt] x [0 ,2 n]. (8) 

We utilize the unitary representation of SU (2) on Fock space 
given by 

R{£2) = e~* aJz e~T>P J ye~T,r J z , ( 9 ) 

where the three components of the angular momentum vector 
J = J( n ) take the second-quantized form 

Ji = Y,(J i '><*P c a c P' ( 10 ) 

ap 

3 For certain symmetries, the compact notation g = \g\e' n is 
schematic as the ’’phase” £2 actually collects several angles. See 
Sec. l2.2l below for the case of SU (2). 

4 Everything discussed throughout the paper can be generalized to 
any (combination of) symmetry(ies) of interest without running into 
any fundamental difficulty. 


with i = x,y,z and (ji) a p = (1 : ce1 7 ,11 : j3). Those one-body 
operators make up the Lie algebra 


•l i ; J j ] L/7, iftJk ! (11) 

k 

where e (; /, denotes the Levi-Civita tensor. The Casimir operator 
of the group built from the infinitesimal generators through a 
non-degenerate invariant bilinear form is the total angular mo¬ 
mentum 

J 2 = E J ?> ( 12 ) 

i=x,y,z 

which is the sum of a one-body and a two-body term, respec¬ 
tively defined as 

A 

jl (\) = E i 2 ( n ) = Ljlp C a C P ’ (13a) 

n=l ap 

A j 

J (2)= E 2'(») •;'(«') = ^ E jjapyScWpC S C r .(13b) 

n/n'=1 ^ aPyS 

Their (direct-product) matrix elements are given by 
Jap — (1 : a l./ 2 |l : P) 

= E <1 : «|yf |1: /3), (14a) 

i=x,y,z 

jjapyS = (1 : «; 2 : / 31 77 11 : y;2 : 8 ) 

= 2 E (l:a|; i |l:y)(2:/3|./,|2:5), (14b) 

i=x,y,z 


from which antisymmetrized matrix elements are obtained through 
JjaPyS = jjapyS ~ jjapdy 

Matrix elements of the irreducible representations (IRREPs) 
of SU (2) are given by the so-called Wigner Z)-functions ll33l 

(^ m \R(£ 2)\^: k ) = 8^8 jj,D j mk (£2), (15) 

where ) is an eigenstate of J 2 and7, 

J 2 \^ M ) = J(J+])h 2 \W 2M ) , (16a) 

Jz\% M ) = Mh\^ M ) . (16b) 

with 2 J £ N, 2 M £Z, J —M £ N and —J < M < +J. The (27+ 

1 (-dimensional IRREPs are labeled by 7 and are spanned by 
the {|1^ M }} for fixed 7 and fi. 

The volume of the group is 

v su( 2 ) = [ d£2 = \6n~, (17) 

J D SU( 2 ) 

such that the orthogonality of Wigner D-functions reads 

Id£2 D J ff K {£2 ) D j m , k , (£2 ) = 8 Jf S MM , 8 kk , . (18) 

^ D su( 2) 27 + 1 

Focusing on the restoration of SU (2) symmetry, it is unnec¬ 
essary to consider the breaking of f/(l) symmetry in the first 
place. We thus limit ourselves for simplicity to (normalized) 
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reference states of the Slater determinant type , i.e. to the case 
where the Bogoliubov reference states defined through Eqs. [2]- 
[3 reduce to 

A 

i^te))=n«« te)+ io)> (i9) 

i= 1 

such that the energy kernel depends solely on the normal den¬ 
sity matrix in this case, i.e. h{g'\g) = h[p 8g ]. The A occupied 
(hole) orbitals in |<J>(,g)} are labeled by (i,j,k,l.. .) while un¬ 
occupied (particle) orbitals are labeled by (a,b,c,d.. .). Greek 
labels (a,/3,y,5 ...) represent occupied or unoccupied states 
indifferently. 


2.3 SR and MR implementations of the EDF method 


The EDF method is embodied in two successive levels of im¬ 
plementation. The first step makes use of the sole diagonal part 
of the kernels. The energy kernel H(g;g) = h(g\g) at play in¬ 
volves one symmetry-breaking state at a time such that this im¬ 
plementation is denoted as the single-reference (SR) level. The 
reference state and the energy are obtained by minimizing the 
latter under the constraint that the magnitude of the order pa¬ 
rameter is fixed to a given value |g|, i.e. the SR energy is given 

by = Min { | € , w>} |^|| along with 


4?i -i [\s\ - l( < £l?)|G ! |‘f>( < ?))| 


( 20 ) 


where G denotes the operator characterizing the order param¬ 
eter, i.e. an operator whose average value is zero in a symmetry- 
conserving state. The minimization of <£j g | leads to solving equa¬ 
tions of motion of the form 

h W (p ( a g) = E l a <Pa ‘’, (21) 


where the diagonal one-body field 


h(*) = 


8p g g* 


( 22 ) 


is obtained as the functional derivative of the diagonal ker¬ 
nel with respect to the diagonal density matrix. The latter is 
characterized by p 8g p = 8 a p, where n\ 8 ^ = 1 for hole states 

and njf I = 0 for particle states. The independence of the energy 
with respect to Q denotes the existence of a pseudo-Goldstone 
mode (whenever g| f 0) and relates to the independence of 
the kernels under a simultaneous rotation of both left and right 
states about the same angle J35). This property indicates that 
off-diagonal kernels must only depend on the difference of the 
angles defining the left and right states, i.e. 

= jvIs'IM^ —f2'), (23a) 

h(\ g ’lQ';\g\,n) = h(\ g '\A\g\,Q-Q l ) 

= h\ g '\\ g \(n-n'). (23b) 


The SR description provides a first account of ground-state 
properties along with a selected set of spectroscopic informa¬ 
tion ED- Allowing the reference state to break symmetries of 


the underlying Hamiltonian is key to incorporate static collec¬ 
tive correlations and thus to address all nuclei, irrespective of 
their closed- or open-shell character. 

The second step makes a full use of the off-diagonal ker¬ 
nels constructed from all pairs of reference states belonging to 
a given set and is thus denoted as the multi-reference (MR) 
level. In state-of-the-art calculations, the set may contain up to 
about 10 9 different symmetry-breaking product states. A key 
feature of the MR description relates to the inclusion of quan¬ 
tum collective fluctuations. In particular, mixing states span¬ 
ning Dc 0 amounts to restoring the symmetry that was possibly 
broken at the SR level. Focusing on such a featui'iQ and stick¬ 
ing to SU(2), the MR energy associated with a state carrying 
good angular momentum reads as 

y _LK'Kf J K U J Kl Dsm dnD J K 1 K (Q)h(Q)N(n) 

EmR = f J K fn su(2) ^ Dp K (Q)N(n) ■ (24) 

To better appreciate the content of Eq.[24] one must realize that, 
as functions defined on D su ^ 2 ), off-diagonal norm and energy 
kernels can be expanded over the IRREPs of the group accord¬ 
ing *0 

N(Q)= £ N J K , K D{ IK (n), (25a) 

JK'K 

h{Q)N(£2) = £ E j k , k N j k , k D j k , k (Q). (25b) 

JK'K 

Exploiting the orthogonality of the IRREPs (Eq.[T8 ]>. the symmetry- 
restored energy can be written as 


F J — 
£ mr — 


I.K'Kf J K * fj^K'K^K'K 

LK'KfPUK'K 


(26) 


and is nothing but a normalized mixing of the coefficients ap¬ 
pearing in expansion[25] In Eqs.[24]and[26] the sum over (K',K) 
mixes the components of the targeted IRREP to remove the 
nonphysical dependence on the orientation of the deformed ref¬ 
erence state. The coefficients of the mixing fj, are generally un¬ 
known and are typically determined utilizing the fact that the 
ground-state energy is a variational minimum. This eventually 
leads to solving a Hill-Wheeler-Griffin equation 1134113011191 


E (E 


K=-J 


J 

K'K 


£ mr) N 


J 

K'K 


f J K = 0 - 


(27) 


Eventually, the MR description not only refines properties al¬ 
ready computed at the SR level but also enlarges the number 
of accessible characteristics of the system, especially by reli¬ 
ably addressing observables that intimately depend on the ful¬ 
fillment of symmetry selection rules, e.g. electromagnetic tran¬ 
sitions of specific multipolarity. All in all, the concept of sym¬ 
metry breaking and restoration is at the heart of the nuclear 

5 The magnitude of the order parameter |g|, i.e. the ’’deformation”, 
is thus omitted throughout the rest of the paper. We refer the reader 
to Ref. |261 for an account of the MR formalism that includes the 
treatment of collective fluctuations associated with |g|. 

6 When treating the restoration of neutron or proton numbers asso¬ 
ciated with the 17(1) group, the corresponding expansion is nothing 
but the Fourier expansion. 
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EDF method as it constitutes a powerful tool to include collec¬ 
tive correlations that are otherwise extremely costly to grasp 
within a symmetry conserving approach. 

2.4 Explicit formulation of the kernels 

Working under the hypothesis that left and right states solely 
differ by a rotation R(Q ), we write 


A 


i*(o)> s fNi°>. 

i=i 

(28a) 

i*(fl))=n4io>. 

i=i 

(28b) 


where rotated orbitals are defined from unrotated ones through 

\cc)=R(Q)\a) = Y,Rp a mP), (29a) 

P 

al=R(Q)al l R\n)=Y J Rpa(nWp, (29b) 

P 

with R a p(Q) = (a\R(Q)\j5) the unitary transformation matrix 
connecting the rotated basis to the unrotated one. The overlap 
between the two Slater determinants |<f>(0)) and |<f>(f2)) can 
be expressed as f45l 

(<f>(0)|<f>(T2)) = detM(Q ), (30) 


kernels that are strictly based on an effective Hamiltonian oper¬ 
ator. While the former route has been the first one followed ll 2 Tl 
I22ll23l and still offers some opportunities l37l . the latter is now 
becoming popular 138II39II40II41II42II43I . 

By construction, off-diagonal energy kernels strictly based 
on an effective Hamilton operator are free from any pathology 
and thus safe to be used meaningfully in MR calculations. Such 
a scheme starts by introducing a tractable effective Hamiltonian 
operator typically containing, e.g., (simple) two and three-body 
pseudo potentials 

H e ff = T + V e ff + W e ff (31a) 

= E f «j3 c a c /3 ( 31b > 

a/3 

+ (ir) E fp r s c aclc s c Y (3lc) 

V z -/ aPyS 

+ fi) E ^y5e? c a44 c ? c e c 5> ( 31d ( 

<xpy8e£ 

where t a p denotes matrix elements of the one-body kinetic en¬ 
ergy operator whereas v^yg and ^^Py8e^ represent antisym¬ 
metrized matrix elements of two- and three-body pseudo poten¬ 
tial operators, respectively. Taking the straight matrix element 
(i.e. mean-field approximation) of the off-diagonal en¬ 
ergy kernel reads, by virtue of the generalized, i.e. off-diagonal, 
Wick’s theorem (GWT) P4l . as 


where M L j (Q ) is the Ax A reduction of R a p (Q ) to the subspace 
of hole states of |<J>( 0 )}. 

As Eq. Q] testifies, the off-diagonal norm kernel is explic¬ 
itly given as the plain overlap between the two reference states 
involved. Formulated phenomenologically, the energy kernel 
is either taken as the off-diagonal matrix element of an effec¬ 
tive Hamilton operatoiQ or as a more general functional of the 
off-diagonal density matrices. In the former case, it takes the 
form of a low-order polynomial functional of the density matri¬ 
ces with specific inter-relations between its coefficients. In the 
latter case the functional typically involves further dependen¬ 
cies on non-integer powers of the off-diagonal density matrices 
and less constrained inter-relations between the coefficients in¬ 
volved. 

In practice, all modern parametrizations characterized by a 
good enough performance (irrespective of the merit function 
used to judge this performance) do not fall in the category of 
functionals formulated as the strict off-diagonal matrix element 
of an effective Hamilton operator. As a result, MR calcula¬ 
tions, although characterized by an apparent success through¬ 
out the first years of their applications, have been shown to be 
plagued with critical pathologies |20l . As of today, the con¬ 
sensus among practitioners is that those pathologies must not 
be overlooked and require a decisive solution ll36l . Two pos¬ 
sible options are in sight; i.e. (i) regularize the pathologies a 
posteriori for a given functional or (ii) limits one-self to energy 

7 In the present context, the wording effective forbids the possibility 
that the Hamiltonian depends on the many-body solution via, e.g., a 
dependence on the (off-diagonal) density (matrix) of the system. 


h(Q) = 



(32a) 

<<P( 0 )|<P(fl)> 

L'aPPpa 

(32b) 

a/3 


5 E 

(32c) 

z a/3 yS 


1 V w eff n 0n On Oil 

^ Zj w aPy8e£ P8a Pep P^y ' 

(3 2d) 


Consequently, h ( Q ) takes the form of a functional of the off- 
diagonal normal density matri^B containing linear, bilinear and 
trilinear term'Q with specific interrelations between the coeffi¬ 
cients of the polynomial. 

Energy kernels of the form given by Eas. [3Tll32l and based 
on extended Skyrme l4Tj and/or finite-range II40II43I pseudo¬ 
potentials (without any density-dependent coupling) are cur¬ 
rently being constructed and implemented. Although it is too 
early to declare success or failure, it already appears that the 
restrictive definition of the energy kernel may result in a lack 
of flexibility I39li36l . As a matter of fact, and as can be ex¬ 
pected from the pioneering work carried out in the 1970’s ED, 
it is a challenge to describe quantitatively all the desired phe- 


s When considering general reference states of the Bogoliubov 
type, additional terms depending on the anomalous off-diagonal den¬ 
sity matrix are generated from the same pseudo-potentials 126 . 

9 This can obviously be extended to a quadrilinear terms by using a 
four-body pseudo-potential etc. 
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nomenoloinPl on the basis of an uncorrelated energy kernel 
(Eq. I32al) . independently of the (tractable) form of //,.nF1. Of 
course, this difficulty must be put in perspective with the in¬ 
valuable capacity to perform safe MR calculations with such 
parametrizations of the energy kernel. If this limitation happens 
to be too significant, one must question the routes that remain 
to be followed. In doing so, one can typically think of 

1. Adding higher-order many-body terms to H e g- to simulate 
missing correlations in the kernel, 

2. Abandoning MR calculations and come back to general 
empirical parametrizations of the diagonal energy kernel. 

It is not clear how to proceed systematically with the first so¬ 
lution that might anyway become quickly impractical. As for 
the second solution, it is not something that one should be will¬ 
ing to follow yet. Indeed, the potentiality of state-of-the-art MR 
codes is too fantasticQ to give up on them at this point in time. 
The goal of the present document is to propose an alternative 
route that combines two wanted features, i.e. an approach that 

1. Follows the strict Hamiltonian-based method and thus leads 
to safe MR calculations by construction, 

2. Provides a way to encode missing correlations into the off- 
diagonal kernel(s) in a controlled and systematic fashion. 

As such, the strategy is to balance the complexity between the 
effective Hamiltonian used on the one hand and the functional 
form of the kernels on the other hand. This is done by building 
off-diagonal norm and energy kernels beyond the uncorrelated 
form given, respectively, by Eas.lTland l32al through a consistent 
expansion of fully correlated off-diagonal many-body kernels. 


3 MBPT of off-diagonal kernels 


the Hamiltonian, // e ff, in terms of which the new EDF kernels 
are expressed. The present discussion is limited to a two-body 
pseudo potential V e g-, i.e. multi-body operators of higher rank 
in /y c ff are discarded for now. Derivations are omitted and the 
interested reader is referred to Ref. Ii28l for technical details. 


3.1 Fully correlated off-diagonal kernels 

We introduce the evolution oncratoiFl in imaginary tirmQ 

(t) = (33) 

and the time-evolved many-body state 

l'I'(T)) ee ^(t) 10(0)), (34) 

which satisfies the time-dependent Schrodinger equation 

H efl \'F(x)) = -d T \'P(x)). (35) 

Having |*P(t)) at hand, we introduce a set of fully correlated, 
off-diagonal and time-dependent kernels 

Af(T,fl) = (*P(T)|l|4>(fl)), (36a) 

H{x,£l) = (W(z)\H eB \0{Q)), (36b) 

Ji(x,n) = (W(T)\Ji\^(n)), (36c) 

J 2 ( x,£2) = ('F(x)\J 2 \0(Q)), (36d) 


which relate to norm, energy, angular momentum projections 
and total angular momentum kernels, respectively. The energy 
and total angular momentum kernels can be further split into 
their one- and two-body components according to 


Building on an earlier work l47l . ab initio MBPT and CC the¬ 
ories based on a symmetry-breaking reference state have been 
generalized in Ref. If28l in such a way that the symmetry is ex¬ 
actly and consistently restored at any truncation order. While 
Ref. (28ll focused on the (breaking and the) restoration of SU (2) 
symmetry, the caseofU(l) symmetry has been formulated even 
more recently (48). This set of novel many-body theories rely 
on the expansion of fully correlated off-diagonal norm and en¬ 
ergy kernels and are, as such, perfectly suited to guide the sys¬ 
tematic construction of their effective EDF counterparts. 

Results obtained in Ref. If28ll are recalled here in the simple 
case of second-order MBPT. In view of their use as a basis for a 
novel EDF scheme, we make explicit the effective character of 

10 This typically includes empirical characteristics of the nuclear 
equation of state, nuclear ground-states observables, including pair¬ 
ing properties, along with gross spectroscopic features of a (large) set 
of nuclei. 

11 In fact, this difficulty is what triggered the use of density depen¬ 
dencies that is now the standard in almost all modern parametrizations 
of the energy kernel. 

12 See e.g. Ref. 03 for the recent development of MR calculations 
of odd nuclei or Ref. m for the computation of isospin-symmetry- 
breaking effects in superallowed Fermi beta decay in view of testing 
the unitarity of the Cabibbo-Kobayashi-Maskawa flavour-mixing ma¬ 
trix. 


H(x,£2) = T(x,£2)+V(x,£2), (37a) 

j 2 (x,n) = j 2 w (x,n)+jf 2 ) (T,n). (37b) 

In the following, a generic operator is denoted as O while its 
off-diagonal kernel is referred to as 

O(x,Q) = {'P(x)\O\0(Q)), (38) 


and is noted in the infinite time limit as 0(°°,£ 2) = 0(£2). Ad¬ 
ditionally, use is made of the reduced kernel defined through 


&(x,£2) 


0{x,£l) 
N(x, 0) ‘ 


(39) 


The advantage of introducing reduced kernels relates to the 
possibility to work with intermediate normalization at rotation 
angle £2 = 0, i.e. with V(x,0 ) = 1 for all x. 

13 Reference f28l being dedicated to ab initio calculations, the for¬ 
malism is formulated having a realistic nuclear Hamiltonian H in 
mind. This does not however prevent one from applying the same 
mathematical formalism in terms of an effective Hamiltonian H e ff. See 
Sec. l4.2l for further comments on the connection between both imple¬ 
mentations. 

14 The time T is real and given in units of MeV -1 . 
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One can show ll28l that the exact lowest energy E ] Q of a 
given IRREP is obtained via 

, Zk'K /£ f J K Jd su( 1 ) Op K ( n )•*?(£*) 

E ° ~ LK'Kfpfl Jd su( 2 ) ^d j * k (Q)^(Q) ’ (40) 

where M ! (L 1 ) and -.Y ( £ 1 ) denote th e fully correlated (reduced) 
many-body kernels. From this point on, the challenge is to 
develop a many-body methodology to expand these kernels 
around the reference state |<J>(0)). How to do this is explained 
at length in Ref. Il28l where both the MBPT and CC expansions 
of off-diagonal kernels are fully worked out. We now report on 
the results obtained from the MBPT expansion at second order 
that are relevant to our purpose. 

3.2 Unperturbed system 

The Hamiltonian is split into a one-body part Hq and a residual 
two-body part Hi 

H eS = H Q +H u (41) 

such that Hq = T + U and Hi = V eff — (/, where 

U^uapckp, (42) 

aft 

is a one-body operator that remains to be specified. Typically, 
we will advocate to take Hq as the the sum of the zero and 
one-body part obtained by normal-ordering H e g with respect to 
the state |0(O)) that minimizes the expectation value of H e g- 
under the possible breaking of SU (2) symmetry, i.e. solving 
deformed Hartree-Fock equations for // e ff one obtains 

( 43 ) 

yS 

Correspondingly,//; is the (symmetry-breaking)normal-ordered 
two-body part of V cff . In this context, one can let 0(0)) break 
SU (2) symmetry spontaneously or force it to do so by adding 
an appropriate Fagrange constraint. This means that the prod¬ 
uct state 10 ( 0 )) is potentially not an eigenstate of J 2 and spans 
several IRREPs of SU (2). 

The operator Hq can be written in diagonal form in terms 
of its (deformed) one-body eigenstates 

p) a \ a a, (44) 


potential is chosen to lie in the energy gap separating particle 
and hole orbitals, i.e. 

e a > p and e, < p . (46) 

Excited eigenstates of Hq are obtained as particle-hole excita¬ 
tions of 10(0)) 

\^m=4aialaj...\^(0)), (47) 

with the eigenenergy 

Hq 0't d»> = (eo + e^;)K;(0)), (48a) 

= e aY e b +...— e,— ej —... . (48b) 

Having |0(O)) at hand, the off-diagonal one-body density 
matrix can be written as m 

p on =f^\t ) Mr j \Q)(j\, (49) 

0=1 

such that it acquires the form 

O on = (l' lh 0\ ( 0 0 \ 

P \ 0 0J + l v R(/2)M- 1 (/2) 0 J 

— P°° + P° Q i (50) 

where t hh is the identity operator on the hole subspace of the 
one-body Hilbert space. The genuinely Q. -dependent part p on , 
which only connects particle kets to hole bras, vanishes for 
Q = 0, i.e. p {){) = 0. The above partitioning of the off-diagonal 
density matrix can be summarized by writing its matrix ele¬ 
ments under the form 



Pa/3 

= n a S a p + (1 — n a)nf 3 ■ 

(51) 

Making particle and hole indices explicit, one obtains equiva¬ 
lently 



Pab= 0, 

(52a) 



pT = o, 

(52b) 



p?f 2 = s ij , 

(52c) 



Paf = Paf ■ 

(5 2d) 


3.3 Second-order expansion of the energy kernel 


where the chemical potential li is introduced for convenience. 
Thus, 10(0)) satisfies 


Expanding the evolution operator in powers of Hi under the 
form 


Hq |0(O)) = £o |0(O)), 

A 

£ o = E(e;-/i). 

i= 1 


(45a) 

(45b) 


The (deformed) Slater determinant |0(O)) necessarily possesses 
a closed-shell character, i.e. there exists a finite energy gap be¬ 
tween the fully occupied shells below the Fermi energy (chem¬ 
ical potential) and the unoccupied levels above. The chemical 


^( T ) = e~ xH ° Te~Jo **i(0 , (53) 

where T denotes the time-ordering operator and where 

Hi (t) = e TH °Hie~ TH ° (54) 

defines the perturbation in the interaction representation, one 
can expand O(Q) in perturbation and represent it diagrammat- 
ically [1281 . This can be achieved by virtue of the GWT and by 
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fjett 
u ct 0~/6 

7 5 




off-diagonal energy kernel, whose diagrammatic representation 
is given in Figs. [3] and |4] is indeed quite rich. To write it in a 
manageable form that is also more amenable to its numerical 
implementation, we first need to express // e ff in convenient left 
and right bi-orthogonal single-particle bases that we now intro¬ 
duce. 

The right basis {| ce)} is obtained by applying the non-unitary 
transformation 


B{Q) = l + p on , (59) 


Fig. 2. First-order Feynman diagrams contributing to li(£ 2). A full 
line carrying an arrow denotes a off-diagonal unperturbed one-body 
propagator G®^(Ti,T 2;T1)- See Ref. ll28l for the details of the dia¬ 
grammatic rules at play. 


onto the original basis {|ce)}. Omitting for notational simplic¬ 
ity the explicit Q dependence of the basis states thus obtained 
and separating original particle and hole states provides 

10 = (60a) 

kc 

|5) = | a ). (60b) 


generalizing the definition of diagonal unperturbed one-body 
propagators to off-diagonal ones 


Particle kets are thus left unchanged. The left basis {(or|} is 
similarly obtained by applying the transformation 


<®|T[ flo (T 1 )aJ(Xi)]|®(fl)) 

- 


(55) 


A key result is the demonstration that the kernel associated with 
any operator factorizes according to 


0(£2) = o(£2)N(£2), (56) 


B~ l (i2) = l-p on , 
onto the original basis {(oc|} such that 

(/I = 01, 

(b\ = (fo|-^(f2 )M w 1 (^)<(|- 

kl 


(61) 

(62a) 

(62b) 


where o(£2 ) denotes the sum of all connected vacuum-to-vacuum 
diagrams linked to the operator O. This writes for the energy 
kernel as 

H(Q) =h(Q)N(Q), (57) 

which thus defines the connected/linked energy kernel h(Q). 

To first ordciFl in MBPT, diagrams contributing to the con¬ 
nected/linked energy kernel are displayed in Fig.[2]and read as 


Hole bras are thus left unchanged. Although we use for sim¬ 
plicity the same notation to characterize states in the left and 
right bases, the tilde is meant to underline their bi-orthogonal 
character. The latter, indicated by («|/3) = 8 a p can be easily 
obtained from B~ l (Q)B(Q) = 1. Given any n-body operator 
O, we introduce the transformed operator O(Q) through 

0{Q)= £ 0 & § {Q)al l ...a^a s ...a r ,(63) 

\ n ' / a..fy...8 


h m (Q) 


(4>(0)\H eS \<P(n)) 

(<p(o)\<P(n)) 

i ia 




2 7 
1 


-eff 
ijc 

-eff ,JH2 


Z ijab 


(58a) 

(58b) 


(58c) 


and is nothing but Eqs. !32bH32cl in which particle and hole 
states have been specified and the off-diagonal density matrix 
has been split according to Eq. [50] The expanded expression of 
the second-order contributions to h(£2), which constitute the 
first corrections to the traditional effective mean-field kernel, 
is too cumbersome to be reported here and we refer the inter¬ 
ested reader to Ref. l28l . This second-order correction to the 

15 The counting of perturbative orders is shifted by one unit relative 
to the convention used in Ref (28). 


where creation and annihilation operators refer to the original 
eigenbasis {|a)} of Ho while left and right indices of the ma¬ 
trix elements refer to the associated bi-orthogonal system in¬ 
troduced above. With these definitions at hand, the first-order 
off-diagonal connected/linked energy kernel (Eq. l58l > can be 
straightforwardly rewritten under the compact form 

h^(Q) = Y,tu(n) + ^vgj(n) (64a) 

i L ij 

= (0(O)|// e ff(f2)|0(O)), (64b) 

which represents the diagonal matrix element of the transformed 
Hamilton operator // e g-(f2) at rotation angle Q. 

The analytic expression of MBPT diagrams making up the 
connected/linked off-diagonal kernel of any operator O can in 
fact be systematically reduced to the expression of the dia¬ 
grams making its connected/linked diagonal kernel, at the sole 
price of involving the transformed operator O(Q) expressed in 
the (T2-dcpcndciit) bi-orthogonal system (28) . Taking second- 
order MBPT as a prime example, and introducing for additional 
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compactness the second-order expression of so-called one- and 
two-body cluster amplitude 




£a ; 


1 -eff 


Y V ijaj Uia 


(65a) 


v ef f 

ijab 


jb 


— e, — ej 


peff 


l]ab ' ' e a + e b -ei-ej 

one obtains the remarkable identit\TH 


P° b f, (65b) 


(65c) 


^{Q) = 


(0(o)| [l + (2) (X2) + £rj { 2 \n) ] o\0{Q)) c 


(0(O)|0(X2)) 

= <0(0)1 [l + ^ { 2 \ q ) + O(fl)|0(O)) c , 


-MeC T\ 



e 



Fig. 3. Second-order Feynman diagrams contributing to t(£ 2). 


where the last matrix element can be worked out on the basis 
of the standard, i.e. diagonal, Wick theorem |49l . Applying this 
to the energy kernel puts us in position to provide its complete 
expression at second order in MBPT under the compact form 

h^(Q) =£%(X2)+£^] (2) (X2)fe(r2) (67a) 

i ia 

ij ija 

(67b) 

i jab 

where Eqs. I67al and I67bl denote the effective kinetic energy 
(Fig. [3]) and potential energy t'( 2 )(T2) (Fig. O parts, re¬ 
spectively. We wish to insist on the fact that, while Ea.[67loffers 
a compact and manageable expression of M 2 )(i2), expanding 
the energy kernel fully would reveal its much richer content 
compared to its first order, i.e. effective mean-field, counter¬ 
part. 


3.4 Consistent expansion of the norm kernel 

So far, we have focused on the expansion of the connected/linked 
part of the kernel associated with an operator O. Dealing with 
the norm kernel amounts to taking 0=1, for which the con¬ 
nected/linked part is trivially equal to one such that one is left 
with the norm kernel itself (see Eq. [56l) by definition. A di¬ 
rect perturbative expansion of N(£ 2) is possible but does not 
constitute the appropriate approach in the present context Il28l . 
The consistent determination of the norm kernel relies on a 
specific argument, i.e. on the fact that applying the symmetry 
restoration scheme to the Casimir (O = J 1 ) of SU{2) and to 
an infinitesimal generator (O = J Z ) that commutes with it must 


Tl 



£ 



Fig. 4. Second-order Feynman diagrams contributing to v(i2). 


16 Choosing U = U HF (Eq. [43} cancels out the first contribution to 

^ fl t(2) (i2) (Eq.m- 

17 The index c indicates the connected character of the matrix ele¬ 
ment. 


give the expected values (respectively J(J + l)7r and MK) in¬ 
dependently of the order at which the expansion is truncated. 
As demonstrated in Ref. (28), this fundamental feature requires 
the reduced norm kernel JV ( Q) to satisfy three coupled ordi- 
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nary differential equations 

-^jV(n) = - l -h{n)jr(n), ( 68 a) 

= +^sin(a) j x (£2 ) jY{£2 ) 

--cos (a)j y {£2)J^{£2), ( 68 b) 

n 

r^-^V(£2) = — ^smtp) cos(a) j x (Q)<yy(Q) 
dy Ti 

- 1 - sin(/3) sin(a) j y {£2)^V{£2) 

Tl 

- l -cos(P)j z (£2)j!£(£2), (68c) 

where { j x {£2 ), j y {£2 ), j z {£2 )} denote the connected/linked ker¬ 
nels of the three infinitesimal generators {J x ,J y ,J z )} of SU (2) 
and whose initial condition is given by the intermediate nor¬ 
malization condition NK(0) = 1. Working at a given MBPT 

order n, one must thus compute {j x '\£2),j y n ^ (£2)J Z '\£2)} 
and solve Eq.[ 68 ]to determine (£2). At second-order, this 
amounts to computing 

jfw = L(AWfl)+L^ t(2) (fl) Ukhm (69) 

i ia 

and solving Eq.[ 68 ]to obtain consistently 

^ 2 \Q) = X {2) {Q)(®(0)\<P{n)). (70) 

In Eq.[70] K W (£2 ) refers to the part of the norm kernel that fac¬ 
torizes in front of the plain overlap. If further limiting Eq. [69] 
to first order, Equation[ 68 ]reduces to the ODEs known l50ll5ll 
to be fulfilled by the plain mean-field overlap, i.e. K'^(X 2 ) = 
1 and yK^^(£2) = (<P(0)l<P(£2)). Eq. [ 68 ] thus constitutes the 
proper generalization, to any order in the many-body expan¬ 
sion, of the ODEs known to be fulfilled by the uncorrelated 
norm kernel. 


3.5 Diagonal kernels 


The many-body scheme developed in Ref. |[28l for off-diagonal 
kernels provides a safe constructive ansatz and avoids the dan¬ 
gers of building off-diagonal kernels as empirical extension of 
known diagonal kernels. Of course, the corresponding many- 
body techniques are compatible with the standard techniques 
applicable to diagonal kernels l45l . i.e. setting £2 = 0 in the 
expansion of off-diagonal kernels provides the same expansion 
for diagonal kernels as standard (symmetry-unrestricted) meth¬ 
ods. The compact form given in Eq. [67] makes quite straight¬ 
forward to see that one does indeed recover standard second- 
order (symmetry-unrestricted) MBPT for £2 = 0. Expanding 
the cluster amplitudes and using U = t/ m for simplicity leads 
to ( 0 ) = 0 for all (a. i) and to 


A P) ( 0 )=J> + 

i 



v e T 

V IJIJ 



K b \ 2 


ea + eb — e, — ej 


, (71) 


along with <jV(0) = 1. This reduction to standard second-order 
MBPT originates from p 00 = 0 or equivalently to the fact that 
the bi-orthogonal system reduces to the mere eigenbasis of // (l 
for £2 = 0, i.e. 7? e ff(0) = // c fi - Although the compact form of 
h£ 1 ' l (£2) connects very naturally with /z^ 2 *( 0 ), it hides the fact 
that its inherent complexity and richness would have forbidden 
to guess its form a priori, to obtain it from h 1 ' 2 ' 1 ( 0 ) via some sort 
of reverse engineering. As such, it was mandatory to develop a 
many-body expansion techniques of genuine off-diagonal ker¬ 
nels from which diagonal ones could be recovered as a partic¬ 
ular case. Eventually, Eq. [Ti] stresses that the diagonal kernel 
at play in SR calculations is explicitly correlated at the level 
of standard second-order MBPT. This expansion can be clearly 
extended to higher orders, at the price of complying with the 
associated increase of computational costs. 


4 Ab initio-driven EDF scheme 

We thus propose to combine the formalism exemplified in the 
previous section with modern EDF techniques. The resulting 
scheme can be summarized as follows, with the corresponding 
algorithm graphically illustrated in Fig.0 


4.1 Algorithm 

A Reference state 

(a) Solve, e.g., symmetry-unrestricted Hartree-Fock equa¬ 
tions in terms of H e g- in the basis of interest to obtain 
the (deformed) reference state |<J>(0)). This amounts to 
using 

h^(0) = h^[p°°], 

as an input diagonal functional. We denote by Nb = Ni, + 
N p the dimension of the one-body Hilbert space, where 
N/, denotes the number of occupied states of |<P(0)) and 
N p the number of unoccupied states. 

(b) Store single-particle energies {e a } and wave-functions 
{<Pa}- 

B Single-reference calculations 

(a) Build the diagonal one-body density matrix p 0(l along 
with the matrix elements of T, \4ffi-4 and any other ob¬ 
servable O in the eigenbasis {cp a } of Hq. 

(b) Compute from it the diagonal energy kernel at the cho¬ 
sen order n 

E SK = h^(0)=h^[p 00 -,{e a }]. 

Proceed similarly for all the other observables O of in¬ 
terest, i.e. compute 

D SR = o W( 0 )= O W[p 00 ;{ ea }]. 

C Multi-reference calculations 

(a) Discretize the intervals of integration over the three Eu¬ 
ler angles £2 = (a, j3, 7 ). 

(b) For each combination of Euler angles 

i. Build the Nb x Nb matrix R a p{£2) = {a\R(£2)\f5) 
and its N), x N/, reduction Mjj{£2 ) to the subspace of 
hole states of |<J>(0)). Compute the inverse M~' (£2). 
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ii. Build the N p x N), rectangular matrix 

p^(Q) = Y j R ai (Q)M ji \Q). 
i= 1 

iii. Build the bi-orthogonal bases according to Eqs. l60l 
and [62] 

iv. Transform the matrix elements of 7’. and any 

other observable O of interest into the bi-orthogonal 
system to generate the matrix elements of t(£2 ), 
V eff (/2), Jk(£2 ) and 0 (£2), respectively. 

v. Compute and store the off-diagonal linked/connected 
kernels at the chosen order n 

h^(Q) = h^[p on -{e a }}, 
o^ n \Q) = o^[p on -{e a }}, 

(n) 

(c) Using j k _ v v ,(I2) for the discretized values of the Eu¬ 
ler angles, along with the initial condition (0) = 1, 
integrate the three coupled ODEs (Eq. l68l > to obtain 

^ n \Q) = NW[p Oi2 ;{ ea }](4>(O)|0(f2)). 

for each combination of the Euler angles. 

(d) Solve the Hill-Wheeler-Griffin equation to obtain the 
weights f j, ('Eci.IT/Ti. 

(e) Calculate the energy of the yrast states through 

j L k'kIP f J K J Dsm dn D£ K (n)hW(£2) ^»> (£2) 

Emr ~ ^k'k/P fi Jd su(2) ^ op K (n) ,/C(«) (q ) 

Proceed similarly to compute other observables O of in¬ 
terest 


4.2 Discussion 

A few further comments are in order regarding the EDF scheme 
we propose. 

1. The proposed scheme is meant to formalize the EDF method 
as used in nuclear physics on a firmer ground than has 
been done in the past. It makes clear that the nuclear EDF 
method invokes fundamentally a triple expansion of the 
many-body solution. First, it involves an ’’horizontal” ex¬ 
pansion in terms of non-orthogonal symmetry-breaking prod¬ 
uct states that pro vides an efficient way to resum so-called 
non-dxnamicaP 1 correlations through the MR mixing. The 
associated off-diagonal kernels N(£ 2) and h(£ 2) are them¬ 
selves ’’vertically” expanded in terms of particle-hole ex¬ 
citations, which corresponds to the treatment of so-called 

18 Care has to be taken for non-scalar operators such that the associ¬ 
ated MR expression must be properly adapted from the one appropri¬ 
ate to the calculation of the energy ED. 

19 Dynamical and non-dynamical refer to the wording at play in 
quantum chemistry. 


Auxiliary diagonal kernel 


h V °] 

o 


Reference state KB 

Symmetry group IB 

^ hU[ P m }\ 

(mm ) 

/ p 00 J 

H 1K1 

\ 

g = {R{tt);£i€D g } 

Generator(s) = J 

IRREP = D J mk (P) 



Single-Reference 


EsR = h < n) [p 00 ;{e Q }] 
0 SR = o(") [p 00 ; {e Q }] 


Off-diagonal kernels at order n in MBPT 


:*Energy fcM (fi) = [p 0!2 ; {e a }] 

t*Observable o''- 71 - 1 (ft) = <y r ‘- : {e Q }] 

:*Generator(s) j< n) (fi) = j< n) [p° n ; {e Q }]-J 

:*Norm Af (n) (fi) = M {n) [p on ; {e a }] 



Multi-Reference 

■Emr = “j 11 1 

c. 

\ Y.UK fftti f Dc dnD{} K (si)hW(n)tf(- >(fi )) 

--- fJ _> nJM 

J K U MR 


Fig. 5. (Color online) Algorithm at play for the ab initio-based EDF 
scheme. 


dynamical correlations. One merit of the many-body for¬ 
malism developed in Ref. l28l is to propose, for the first 
time, a consistent simultaneous expansion of both dynam¬ 
ical and non-dynamical correlations, including their inter¬ 
ference. Indeed, those two types of correlations are not or¬ 
thogonal, i.e. not mutually exclusive, such that one may 
choose to convert one into the other. The dual expansion 
that lies at the core of the ab-initio-driven MR formalism 
is thus meant to optimize consistently the sharing between 
both kinds of correlations. Eventually, the scheme is imple¬ 
mented on the basis an operator // e ff that must be seen as 
an effective generator of optimal EDF kernels at each given 
order of this dual expansion, i.e. its free parameters are to 
be optimized to reproduce experimental data at that chosen 
order. Whenever augmenting the explicit content of the ker¬ 
nels or/and of // e ff itself, the parameters of the latter must 
be re-optimized fully. At this point, the form of // e ff is to 
be guessed empirically. However, one can hope to design in 
the future a truly systematic low-energy effective theory that 
formulates the consistent expansion of the effective opera¬ 
tor // e ff in connection with the dual expansion of the many- 
body propagation to eventually achieve the triple expansion 
we alluded to above. In doing so, renormalizability of the 
kernels should guide the possible form of // e ff at any given 
order in the many-body expansion, e.g. see Ref. |[55l for 
the discussion dedicated to the second-order diagonal ker- 
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nel in isospin-symmetric infinite nuclcaiFl In such a low- 
energy effective theory, the generator // L ,n would formally 
converge to a realistic nuclear Hamiltonian H at play in ab 
initio calculations whenever the many-body propagation is 
implemented to a sufficiently high order. 

2. The present scheme clarifies that building many-body cor¬ 
rections into the off-diagonal connected/linked energy ker¬ 
nel must be accompanied by a consistent enrichment of the 
off-diagonal norm kernel N(Q). This is at variance with 
usual practice that consists of keeping the norm kernel to be 
the plain overlap between the two reference states involved 
independently of the many-body content of the energy ker¬ 
nel. This crucial point demonstrates that any well-founded 
formulation of the nuclear EDF method must consider en¬ 
ergy and norm kernels on the same footing, which further 
stresses the necessity to focus on off-diagonal kernels. In¬ 
deed, the fact that the diagonal norm kernel can always 
be normalized to 1 , i.e. one can always redefine .A' ( Q) 
in place of the original N(£ 2), has undermined the neces¬ 
sity to understand how the norm kernel behaves as soon as 
the two reference states involved differ from one another. 
As for SU (2) symmetry, the consistency between both ker¬ 
nels is ensured via the ODEs fulfilled by N(£ 2) that involve 
off-diagonal kernels of the infinitesimal generators of the 
symmetry group computed at the same order as the energy 
kernel. For the one-dimensional 1/(1) group, the norm ker¬ 
nel is related to the particle number operator kernel via one 
first-order ODE | |48| . 

3. Mean-field expressions denote very restricted mathematical 
forms of the kernels, which happen to be functionals of the 
sole off-diagonal one-body density matrix, i.e. 

h (l \L 2 ) = h w [p on ], 
o^\£2)=o^[p on ], 

jP(n) = jW[p on ], 

= (&(o)\<P(n)). 

As soon as one includes explicit corrections to them, the 
kernels naturally become functionals of both the off-diagonal 
density matrix and the A), single-particle energies associ¬ 
ated with the unperturbed HamiltoniarFI 

h^(Q)=h^[p° n -,{e a }], 

o^(Q)=o^[p on -{e a }}, 

(n)=jf [p on ;{e a }], 

^(n) {Q ) = ^)[p oii -{e a }}(<P(0)\<P(Q)). 

Consequently, our proposal naturally extends to off-diagonal 
kernels the notion of orbital- and energy-dependent func¬ 
tionals at play in Density Functional Theory applicable to 
electronic systems urn 

20 Renormalizability implies in particular that the results are inde¬ 
pendent of the particular regularization scheme used. 

21 Note that the single-particle energies are invariant under the ro¬ 
tation of Hq, i.e. under the transformation that leads from | <£> (0) ) to 
|<J>(f2)). As such, they equally refer to |<J>(0)) and |<J>(f2)). 


4. Off-diagonal kernels are systematically improvable beyond 
traditional mean-field matrix elements, i.e. they are based 
on a systematic expansion. We advocate here to employ the 
simplest expansion based on low-order perturbation theory 
but other more involved choices could be envisioned, such 
as a low-order coupled-cluster truncation. Eventually, one 
aims at an optimal compromise between the complexity of 
the many-body formalism and that of the effective Hamil¬ 
tonian, enabling the resulting method to be at the same time 
sufficiently rich and computationally manageable over the 
nuclear chart. 

5. The reference state |<J>(0)) is typically to be obtained on the 
basis of a mean-field-like kernel, e.g. the first-order diago¬ 
nal kernel h 11 1 (0). Once this is done, the actual energy of 
the system can be implemented at a higher order «, both 
at the SR and MR levels. While a link is kept between 
both steps in the sense that both /z^^O) and h^ n \Q) are 
generated from the same effective Hamiltonian // e fr, this 
demonstrates that the determination of the reference state is 
(can be) decoupled from the actual computation of the to¬ 
tal energy. This is nothing but the logic typically followed 
in MBPT. As a matter of fact, one could go even further 
and choose to optimize the reference state on the basis of 
a kernel that is entirely disconnected from h^ n \Q.), e.g. 
one could optimize | </>(())) from an auxiliary, e.g. Skyrme, 
parametrization of the diagonal kernel (see step (E) in Fig. [5} 
that effectively account for some correlations before pro¬ 
ceeding to the spuriosity-free and explicitly-correlated cal¬ 
culation of £sr and on the basis of h' n HQ). Such a 
scheme would require that the parameters of // e ff are ad¬ 
justed relative to the auxiliary parametrization used to opti¬ 
mize | <£(0)) such that the parameter space effectively defin¬ 
ing the approach would be larger than in the approach we 
presently advocate to follow. 

6 . Only at lowest order (n = 1) can the symmetry-restored en¬ 
ergy be factorized as the expectation value of // c n in a 
symmetry-projected wave-function l28l . As soon as higher 
orders are included (n > 2 ), this becomes impossible, i.e. 
the underlying symmetry-conserving wave-function is only 
implicit. This is actually necessary to obtain a connected/linked 
expression of the energy kernel, which is itself mandatory 
for the method to be size extensive, i.e. for the energy to 
scale correctly with particle number l54l . 

7. The present scheme offers a way to include many-body cor¬ 
relations consistently in the computation of all the observ¬ 
ables of interest, i.e. not only for the energy but also for 
any observable O such as, e.g., charge radii, electromag¬ 
netic moments etc. 

8 . Deriving the kernels from a true operator via explicit many- 
body techniques offers the possibility to carry out spuriousity- 
free, i.e. safe, MR-EDF calculations. It however does not 
guarantee it, i.e. an incautious truncation of the full many- 
body expansion may induce spurious self-interaction and 
self-pairing contributions I22U26I . In the present discussion 
based on particle-number conserving reference states and 
properly antisymmetrized matrix elements, the procedure 
is always safe. When expanding more general off-diagonal 
kernels defined from Bogoliubov reference states ||48l , more 
attention must however be paid to this question. Eventually, 
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truncation schemes based on strict perturbation theory hap¬ 
pen to be indeed safe. 

9. To reduce the cost of the symmetry restoration (56], it will 
be of interest to apply the topological gaussian overlap ap¬ 
proximation 1571581 to the second-order off-diagonal ker¬ 
nels. 


5 Conclusions and perspectives 

This programmatic document lays down the possibility to build 
novel parametrizations of the off-diagonal energy and norm 
kernels that lie at the heart of the nuclear energy density func¬ 
tional method. The proposal is to exploit the ab initio many- 
body formalism recently proposed in Ref. 1281 to guide the 
construction of safe, explicitly correlated and systematically 
improvable parametrizations. The many-body formalism of in¬ 
terest relies on the concepts of symmetry breaking and restora¬ 
tion that have made the fortune of the nuclear EDF method and 
is, as such, amenable to this guidance. We have detailed the 
proposal in its basic form, including the general equations and 
formulas that need to be implemented. In the mid-term future, 
we wish to investigate this scheme according to the following 
plan 

1. Although St/(2) and Slater determinants have been used 
in the present document as a prime example to expose the 
new EDF scheme, we plan to implement it first in connec¬ 
tion with the breaking and the restoration of t/(l) (particle 
number) symmetry, which requires the handling of Bogoli- 
ubov reference states. The ab initio many-body method that 
can serve as a guidance and that will allow us to fully map 
the present discussion to the U( 1) group has been worked 
out in Ref. l48l . 

2. Although this is a long-term objective, we will not exploit 
renormalizability as a guiding principle to build H e g at first. 
We will use two-body contact interactions containing both 
central and spin-orbit parts and regulated either with a Gaus¬ 
sian form factor or a sharp cut off. 

3. The kernels will be computed at second order in perturba¬ 
tion theory. Although the computation of n = 2 kernels is 
more costly than that of traditional mean-field (n = 1) ker¬ 
nels, it is the least costly many-body method to correct for 
the deficiencies of the latter. As a matter of fact, there ex¬ 
ist methods developed in quantum chemistry to reduce the 
corresponding cost from a scaling that naturally goes as Nl 
down to a scaling that goes as, e.g., N? x N/, ll59l . In consid¬ 
ering these (approximation) techniques, it will be essential 
to only retain those that do not induce any pathology as¬ 
sociated with spurious self-interaction and self-pairing pro¬ 
cesses l22l|26| , which are formally avoided in the first place 
in the proposed scheme. 

The present discussion has focused on the restoration of 
symmetries, i.e. on multi-reference calculations that carry the 
mixing over the angle(s) Q of the order parameter g of the 
broken symmetry. This constitutes only one side of the coin 
of MR-EDF calculations that can/should also treat fluctuations 
over the norm |g| of the order parameter. The possibility to 
formulate the mixing over |g| on the same footing as what 


is proposed here for the symmetry restoration remains to be 
formulated. Although more challenging in several respect, the 
corresponding many-body formalism is currently being formu¬ 
lated (60l and will thus allow a consistent finalization of the 
presently proposed EDF scheme. Once this is done, it will be of 
interest to derive the corresponding (extended) random phase 
approximation. This can be achieved by taking the limit of the 
second-order off-diagonal kernels where \tp( g ■') and dif¬ 

fer harmonically from a common reference state 116 1 Il62ll . 
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